Polaron in t-J model 
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We present numeric results for ground state and angle resolved photoemission spectra (ARPES) 
for single hole in t-J model coupled to optical phonons. The systematic-error free diagrammatic 
Monte Carlo is employed where the Feynman graphs for the Matsubara Green function in imaginary 
time are summed up completely with respect to phonons variables, while magnetic variables are 
subjected to non-crossing approximation. We obtain that at electron-phonon coupling constants 
relevant for high Tc cuprates the polaron undergoes self-trapping crossover to strong coupling limit 
and theoretical ARPES demonstrate features observed in experiment: a broad peak in the bottom 
of the spectra has momentum dependence which coincides with that of hole in pure t-J model. 
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In the context of broad interest to the challenging prop- 
erties of the high-Tc superconductors, a single hole in the 
Mott insulator has been studied extensively in terms of 
the t-J model 
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where cja is projected (to avoid double occupancy) 
fermion annihilation operator, < 2 is the occupa- 
tion number, Si is spin-1/2 operator, and {ij) denotes 
nearest-neighbour sites in two-dimensional lattice. Dif- 
ferent theoretical approaches (for review and recent stud- 
ies see Refs. 2i> 3e> S) gi'^s consistent results for the 
Lehman Spectral Function (LSF) Lk(cj) = —n^^'^Gii{uj) 
of the Green function Gu{uj). Namely, the LSFs at all 
momenta have a quasiparticle (QP) peak in the low en- 
ergy part together with a broad incoherent continuum 
extending up to the energy scale of the order of t. Sharp 
QP peak in the LSF of the ground state at momentum 
k = (7r/2, 7r/2) has the weight Z ^ J/t. This QP peak is 
sharp at all momenta and it's energy dispersion has the 



bandwidth W 



J/t 



J. More advanced t-t'-t"-,l model 



takes into account long range hopping amplitudes t' and 
t" which alter the bandwidth and dispersion of the QP 
resonance [1 B Sll H IHIi^ • K was shown ^ that, de- 
pending on the parameters of the t-t'-t"-J model, the QP 
peak can be either enhanced or completely suppressed at 
a large part of the Brillouin zone. However, at parame- 
ters which are needed for description of ARPES measure- 
ments on carefully studied [13 Sr2Cu02Cl2 {J/t w 0.4, 
t' /t « —0.3, see [aHIlllithe QP peak remains well de- 
fined for all momenta Possible damping of the 

peak 0, in contrast to experiment, is much less than the 
QP bandwidth. Therefore, the properties of QP peak in 
t-t'-t"-3 model at realistic parameters are the same as in 
generic t-J model. 

Experimentally, ARPES in the undoped cuprates re- 
vealed the LSF of a single hole 0, 0| and observed 
dispersion of the lowest peak in the LSF is in good 
agreement with the theoretical predictions of the t-t'-t"- 
J model H II i, H ES im ■ The puzzling point is that. 



in contrast to the theory, experiments never show sharp 
QP resonance and a broad peak with the width of the 
order of 0.1-0.5eV (w t) is observed instead. Note, that 
broadening is seen just in undoped systems where the 
ground state of the single hole is the lowest energy state 
in the Hilbert space of the (N-l)-electron problem. One 
can rule out a possibility of an extrinsic origin of this 
width since the doping introduces further disorder, while 
a sharper peak is observed in overdoped region. 

The role of electron-phonon (e-ph) coupling has gained 
the recent renewed interest. One reason is that the 
ARPES data in doped metallic cuprates revealed the en- 
ergy dispersion strongly renormalized by e-ph interaction 
^1]. Besides, the strong e-ph interaction is crucial for ex- 
planation of renormalization and lineshapes of phonons 
observed in neutron scattering experiments 0, . In 
addition, the large isotope effect on Tc for underdoped 
cuprates and on the superfluid density at th e op timal 
doping suggests the vital role of e-ph coupling . 

Strong and intermediate coupling regimes of e-ph inter- 
action in t-J model can not be studied neither by exact di- 
agonalization ^3 on small clusters nor by self consistent 
Born approximation (SCBA) where both magnon and 
phonon propagators are subjected to non crossing ap- 
proximation (NCA) which neglects the vertex corrections 
E, 113 ■ Small system size implies a discrete spec- 
trum and, hence, neither crossover from large to small 
size polaron nor the problem of linewidth in ARPES 
can be addressed in the former approach. The latter 
method omits the Feynman diagrams with mutual cross- 
ing of phonon propagators and, hence, neglects the vor- 
tex corrections which are crucial for treating the e-ph in- 
teraction in the strong coupling regime (SCR). One can 
use the NCA for the interaction of the hole with mag- 
netic system because spin S=l/2 can not flip more than 
one time and magnon cloud around the hole saturates. 
To the contrary, phonon-phonon non-crossing approxi- 
mation (PPNCA) fails for SCR since number of phonon 
quanta around the hole is not limited. 

The key role of the vortex corrections for e-ph inter- 
action in SCR can be demonstrated by numerically ex- 
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act Diagrammatic Monte Carlo (DMC) method [2l|, 122, 
01 where Feynman graphs for the Matsubara Green 
function of a hole in phonon bath are generated by Monte 
Carlo and summcd up without systematic errors. In the 
framework of this technique one can compare results of 
complete summation of Feynman expansion with those of 
restricted summation where diagrams with inter-crossing 
of phonon propagators are neglected. In order to extract 
the effect of vortex corrections on results we considered 
conventional two dimensional Holstein model, where the 
complete and PPNCA summation of Feynman expansion 
by DMC can be performed exactly. In this model hole 
freely hops with the amplitude t (t is set to unity below) 
and interacts with dispersionless (frequency = const) 
optical phonons by short range coupling 7 



j6k + h.c. 
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Exact and PPNCA results for energy and Z-factor of po- 
laron are in good agreement at small values g < 0.2 of di- 
mensionless interaction constant g = j'^/{8tQ) while are 
crucially different in the intermediate and strong coupling 
regimes. E. g., for fl/t = 0.1 exact result shows crossover 
to SCR g> g'jjK. 0.6 while in PPNCA the polaron is 
in weak coupling regime even for g = 60. We conclude 
that DMC is the only method which can treat interme- 
diate and strong coupling regime of the e-ph interaction 
for the problem of one hole in the macroscopic system. 

In this Letter we present a study of a single hole in t- 
J model interacting with dispersionless optical phonons 
© by DMC IH m lis method. We found that due 
to slowing down of the hole by spin flip cloud around 
it, the hole in t-J model is subject to the stronger e- 
ph coupling than the freely propagating hole and hence 
undergoes the crossover to SCR at smaller coupling. This 
is in contrast to the naive expectation that the small Z 
factor reduces the e-ph coupling in t-J model. Besides, 
we found that SCR occurs at e-ph couplings which are 
typical for high Tc materials. Finally, our results for SCR 
qualitatively reproduce data of ARPES experiments: a 
hroad peak, whose energy dispersion is similar to that of 
pure t-J model, dominates in low energy part of LSF. 

In the Standard spin-wave approximation in momen- 
tum representation (25i | the dispersionless hole eo = const 
(annihilation operator is ft-k) propagates in the magnon 
(annihilation operator is ak) bath 
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-sign{Jk)^/JT^^v^y/J2l^, and vy, = a/1 - 7k- 
We chose the value J/t — 0.3 at which NCA for rnagnons 
is shown to be sufRciently good approximation '25] and 
set the phonon frequency Q — 0.1 to make it similar to 
experimental value. 

We generate Feynman expansion of Matsubara Green 
function of a hole in momentum representation for the 
infinite system at zero temperature by DMC 12^ |^ . 
Then we obtain LSF by stochastic optimization method 
|2^ l26j which resolves equally well both sharp and 
broad features of the spectra. Diagrams with crossing 
of a magnon propagator by both magnon and phonon 
lines are neglected |23 though the diagrams with inter- 
crossing of phonon propagators are taken into account. 
The DMC algorithm for pure t-J model eqs. H3I4() is 
equivalent to the macroscopic limit of SCBA approach on 
finite lattices "25'|. When PPNCA is introduccd to DMC 
at finite couplings to phonons, the algorithm is nothing 
but the thermodynamic limit of SCBA approach on finite 
lattices to the complete model (|2l4f) when both magnons 
and phonons are taken into account in NCA 1201 . Com- 
parison of our data with results of Refs. pol l25j shows 
that Z-factor and energy of the lowest peak are weakly 
influenced by finite size corrections sincc the relative dis- 
crepancy is less than 10~^. The shapes of LSF in SCBA 
approach are very similar to those obtained by DMC. We 
observed only slight discrepancy in widths and energies 
of high energy peaks. Therefore, the main advantage of 
DMC over other existing methods is the possibility to 
take into account phonon-phonon vertex corrections in 
macroscopic system. 
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1: Hole LSF in ground state at J/t = 0.3: (a) for 
0; (b-d) low energy part for g=0.1445 [7 = 0.34] (b), 



g=0.2 [7 = 0.2] (c), and g=0.231125 [7 = 0.43] (d). 



with magnon dispersion Wk = 2J-y/l — 7^, where 7k 
(cos kx + cos ky)/2. The hole is scattered by magnons 
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with the scattering vortex Afk,q — 4t(uq7k-q 
?^q7k), where Uk = a/(1 + i^k)/{2i^k), Uk 



Figure ^ presents our results for LSF in the ground 
state. Three low energy peaks for g = O (Fig. ^) are 
string resonances 12^ since their energies, as we found 
(4) in calculations of LSF for different exchange constants 
0.1 < J < 0.4, obey the scaling low a + bi{J/tf/^ 
(a ^ -3.13, 61,2,3 = 4.89/2.09/1.63). With increasing 
of g these peaks are observed up to 5 = 0.21 (Figs. 
and^) while the peak with highest energy broadens and 
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FIG. 2: Dependence on coupling strength at J/t = 0.3: (a) 
energies of lowest LSF resonances; (b [c]) energy [Z-factor] of 
lowest peak in DMC (circles) and PPNCA (triangles). Lines 
are to guide an eye. The errorbar, if not shown, is less than 
the point size. 
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FIG. 3: LSF of the hole at J /t = 0.3 and g = 0.231125: 
(a) full energy range for k = (7r/2,7r/2); (b-d) low energy 
part for different momenta. Slanted arrows show broad peaks 
which can be interpreted in ARPES spectra as coherent (C) 
and incoherent (I) part. Vertical arrows indicate position of 
ground state resonance which is not seen in the vertical scale 
of the figure. 



disappears at higher couplings (Fig. ^). 

Dependence of the peak energies (Fig.[2^) and ground 
State Z-factor (Fig.|3;) on g resembles picture inherent in 
self-trapping phenomenon [26L [29l | : the states cross and 
hybridize at critical coupling g^ j 0.19 and Z-factor of 
ground state resonance rapidly decreases. We note, that 
PPNCA result does not show transition into SCR even at 
considerably larger g (see Fig. |2l3 and Fig. Accord- 
ing to general understanding of self-trapping crossover, 
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FIG. 4: Dispersion of resonances energies at J/t = 0.3: broad 
resonance (filled circles) and lowest polaron resonance (filled 
squares) sd, g = 0.231125; third broad resonance (open circles) 
and lowest polaron resonance (open squares) at <? = 0.2. The 
solid curves are dispersions of a hole in pure t-J model at 
J/t = 0.3 {Wj/t=o.3 = 0.6): Ermn = -2.396 (e„,„ = -2.52) 
for solid (dotted) line. 



at small g the ground state is weakly coupled to phonons 
while excited resonances have strong lattice deforniation. 
At critical coupling the crossing and hybridization of 
these states occurs and for higher g the roles of these 
states exchange: the lowest state is strongly trapped by 
lattice deforniation while the upper ones are nearly free. 
Therefore, above the crossover point g > g^ j one expects 
that the lowest state should be dispersionless while the 
upper resonances have to show considerable momentum 
dependence. This assumption is supported by momen- 
tum dependence of LSF well above the crossover point 
(Fig. O. The energy of the lowest peak is momentum 
independent while the bandwidth of the upper broad res- 
onance (Fig. ^ is the same (VFj/t=o.3 = 0.6) as it is in 
the pure t-J model. The most surprising peculiarity of 
the momentum dependence of broad resonance is that it 
is exactly the same as expected for t-J model with no 
coupling to phonons: it obeys the scaling relation 



Ek = emm + W j/t{[cOSkx +COS kyf /5 + 

[cos(fcj; + ky) + cos{kx - fcj^)] V4}, 



(5) 



which describes dispersion for pure t-J model in the broad 
range of exchange constants 0.1 < J/t > 0.9 js^]. We 
emphasize, that unrenormalized dispersion of the upper 
resonances is the general property of SCR. E. g., for the 
coupling g = 0.2, which is slightly higher than the critical 
couphng (7j_j, (see Fig.^ and|21), dispersion of the upper 
resonance also obeys eqs. (01 while the lowest peak is 
momentum independent (see Fig.Q. 

The behavior of LSF in SCR is exactly the same as 
it is observed in ARPES experiments. The LSF con- 
sists of broad QP peak and high-energy incoherent con- 
tinuum (see Fig.|31). Besides, momentum dependence of 
the broad QP peak is similar to dispersion of sharp res- 
onance in the pure t-J model (see Fig. 0J. The lowest 
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lying dispersionless peak in SCR has small Z-factor and 
can not be discerned in ARPES experiments: it's spectral 
weight is small (Z ~ 0.015) a,t g — 0.2 (Fig.^) and com- 
pletely suppressed {Z < 10'^) at g ^ 0.231125 (Fig.Uli). 
On the other hand, the momentum dependence of spec- 
tral weight Z' of broad dispersive resonance in SCR is 
akin to that in pure t-J model. For g = 0.231125 the 
weights of broad peak at k = (7r/2,7r/2) 0.27) and at 
k = (0, 0) (~ 0.05) coincide with Z-factors of sharp reso- 
nances at corresponding momenta in pure t-J model. Our 
calculation show that this similarity is the robust feature 
of SCR: ^(^/2,^2) ~ 0.27±0.01 and Z[o ~ 0.05±0.005 
in the broad range of coupling constants 0.21 < g < 0.27. 

Comparing the critical interaction g^ j « 0.19 for the 
hole in t-J model and critical coupling gj^ « 0.6 for Hol- 
stein model with the same value of hopping t, we conclude 
that interaction with spins enhances e-ph coupling and 
accelerates transition into SCR. Emission of magnons on 
each hopping shrinks the bandwidth of lowest lying res- 
onance by the factor of J/t while the reduction of e-ph 
vertex by the factor Z seems to be absent. This makes 
the influence of e-ph interaction on low energy part of 
LSF more effective. We emphasize that the critical inter- 
action for transition to SCR is low enough for t-J model 



to bring the realistic system into SCR. E. g., theoretical 
estimate of the interaction strength from the fitting of 
phonon energies and lineshapes to neutron scattering cx- 
periment gives rather large magnitude for e-ph coupling 
m La(Sr)Cu04 El El. Since the interaction vertex is 
momentum dependent in the model (Tsl ITl , we can es- 
tablish only lower and upper bound for effective coupling 
g in our model with short-range interaction. The av- 
eraging over the Brillouin zone gives the lower bound 
5rap — 0-15 while the maximal value at the Brillouin 
zone boundary determines the upper bound .g™p'^ < 0.5. 
However, since self-trapping is governed mainly by short- 
range coupling pgj . effective constant in our model is 
more close to the upper bound g^^- Thus, realistic high 
Tc cuprates are expected to be in SCR. 

Finally, we conclude that puzzling behavior of ARPES 
spectra in undoped high Tc materials, which manifests 
oneself in unexpectedly broad quasiparticle peak with dis- 
persion corresponding to pure Mott insulator model, is 
driven by strong electron-phonon interaction. 
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